Negative continuum effects on the two— photon decay rates of hydrogen— like ions 
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Two-photon decay of hydrogen-like ions is studied within the framework of second-order per- 
turbation theory, based on relativistic Dirac's equation. Special attention is paid to the effects 
arising from the summation over the negative-energy (intermediate virtual) states that occurs in 
such a framework. In order to investigate the role of these states, detailed calculations have been 
carried out for the 2si/2 — > lsi/2 and 2pi/2 —* lsi/2 transitions in neutral hydrogen H as well as 
for hydrogen-like xenon Xe^^+ and uranium U ions. We found that for a correct evaluation of 
the total and energy-differential decay rates, summation over the negative-energy part of Dirac's 
spectrum should be properly taken into account both for high-Z and \ow-Z atomic systems. 
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I. INTRODUCTION 



Experimental and theoretical studies on the two- 
photon transitions in atomic systems have a long tra- 
dition. Following seminal works by Goppert-Mayer [l| 
and by Breit and Teller [2| a large number of investiga- 
tions have been performed in the past which focused on 
the decay of metastable states of light neutral atoms and 
\ow-Z ions. These investigations have dealt not only with 
the total and energy-differential decay rates fl, 0,01 but 
also with the angular distributions p, 0, B U and even 
polarization correlations between the two emitted pho- 
tons [1^, [lllj ■ Detailed analysis of these two-photon 
properties have revealed unique information about elec- 
tron densities in astrophysical plasmas and thermal X- 
ray sources, highly precise values of physical constants 
[13| , structural properties of few-electron systems includ- 
ing subtle quantum electrodynamical (QED) effects 
as well as about the basic concepts of quantum physics 
such as, e.g., non-locality and non-separability 

Beside the decay of metastable states of low-Z sys- 
tems, much of today's interest is focused also on the 
two-photon transitions in high-Z ions and atoms which 
provide a sensitive tool for improving our understanding 
of the electron-photon interactions in the presence of ex- 
tremely strong electromagnetic fields [l6| . In such strong 
fields produced by heavy nuclei, relativistic and retar- 
dation effects become of paramount importance and may 
strongly affect the properties of two-photon emission. To 
explore these effects, therefore, theoretical investigations 
based on Dirac's equation have been carried for the total 
and energy-differential decay rates [13, [H, [H, [13, [Mj 
as well as for the angular and polarization correlations 
[13, [H, I13|- In general, relativistic predictions for the 
two-photon total and differential properties have been 



found in a good agreement with experimental data ob- 
tained for the decay of inner-shell vacancies of heavy 
neutral atom s [25l . [2q and excited states of high-Z few- 
electron ions [231 . 

Although intensive experimental and theoretical efforts 
have been undertaken recently to understand relativistic 
effects on the two-photon transitions in heavy ions and 
atoms, a number of questions still remain open. One of 
the questions, which currently attracts much of interest, 
concerns the role of negative energy solutions of Dirac's 
equation in relativistic two-photon calculations. Usually, 
these calculations are performed within the framework of 
the second-order perturbation theory and, hence, require 
summation over the (virtual) intermediate ion states. 
Such a summation, running over the complete spectrum, 
should obviously include not only positive- (discrete and 
continuum) but also negative-eigenenergy Dirac's states. 
One might expect, however, that since the energy release 
in two-photon bound-bound transitions is less than the 
energy required for the electron-positron pair produc- 
tion, the contribution from the negative part of Dirac's 
spectrum should be negligible even for the decay of heav- 
iest elements. From practical viewpoint, this assumption 
justifies the restriction of the intermediate-state summa- 
tion to the positive-energy solutions only. Exclusion of 
the negative continuum would lead, in turn, to a sig- 
nificant simplification of the the second-order relativis- 
tic calculations especially for many-electron systems for 
which the problem of (many particle) negative continuum 
still remains unsolved. 

Despite the (relatively) small energy of two-photon 
transitions, the influence of Dirac's negative continuum 
in second-order calculations should be further questioned 
because of possibility for production and subsequent an- 
nihilation of the virtual anti-particles. It has been ar- 
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gued, for example, that transitions involving positron 
states have to be taken into account for the proper de- 
scription of Thomson scattering [2^ , interaction of ions 
with intense electromagnetic pulses |29l . Isoj in the "un- 
dercritical" regime as well as magnetic transitions in two- 
electron ions [31I [H, [3^. Moreover, the first step to- 
wards the analysis of negative-energy contributions to 
the two-photon prop erties has been done by Labzowsky 
and co-workers [3J| who focused on ElMl and E1E2 
2pi/2 ~^ lsi/2 total decay probabilities. The relativistic 
calculations have indicated the importance of negative- 
energy contributions not only for high-Z but also for 
low-Z hydrogen-like ions. 

In this work, we apply the second-order perturbation 
theory based on relativistic Dirac's equation in order to 
re-analyze atomic two-photon decay. We pay special 
attention to the influence of negative continuum solu- 
tions on the evaluation of the transition amplitudes and, 
hence, on the total and energy-differential decay rates. 
For the sake of clarity, we restrict our analysis to the 
decay of hydrogen-like ions for which both the positive- 
and negative-energy parts of Dirac's spectrum can be 
still studied in a systematic way by making use of a finite 
basis set method [l^. Implementation of this method 
for computing relativistic second-order transition am- 
plitudes is briefly discussed in Sections III Al and III Bl 
Later, in Section III CI we consider an alternative, semi- 
classical, approach which allows analytical evaluation of 
the negative-energy contributions to the two-photon ma- 
trix elements and transition rates. These two — semi- 
classical and fully relativistic — approaches are used in 
Section HIT] to calculate the energy-differential and total 
decay rates for several multipole terms in the 2sx/2 ~* 
Is 1/2 and 2pi/2 — > Is 1/2 two-photon decay of neutral 
hydrogen as well as hydrogen-like xenon Xe^^^ and ura- 
nium U^^"'" ions. Based on the results of our calculations, 
we argue that both the total transition probabilities and 
the photon energy distributions can be strongly affected 
by the negative-state contributions; this effect is most 
clearly observed for the non-dipole transitions not only 
in high-Z but also in (non-relativistic) low-Z domain. 
Brief summary of these findings and outlooks are given 
finally in Section [TVl 



II. THEORY 

A. Differential and total decay rates 

Not much has to be said about the basic formalism 
for studying the two-photon transitions in hydrogen-like 
ions. In the past, this formalism has been widely applied 
in order to inve stig ate not only the total decay proba- 
bilities [13, m, [ii) [13 but also the energy as well as 
angular distributions [23| and even the correlation in the 
polarization state of the photons [l^ [2^ . Below, there- 
fore, we restrict ourselves to a rather brief account of the 
basic expressions, just enough for discussing the role of 



negative-energy solutions of Dirac's equation in comput- 
ing of the two-photon (total and differential) rates. 

The properties of the two-photon atomic transitions 
are evaluated, usually, within the framework of the 
second-order perturbation theory. When based on 
Dirac's equation, this theory gives the following expres- 
sion for the differential in energy decay rate: 



dw 
dwi 



_UJiUJ2_ 

(27r)3c3 



E 



{f\A*2\i^) {i^\At\i) 



E,, — Ej 



{f\Al\,y) {iy\A*\z) 



E„ — Ei 



- UJ2 



(1) 



where the transition operators A* with j — 1,2 describe 
the (relativistic) electron-photon interaction. For the 
emission of photons with wave vectors kj and polariza- 
tion vectors ej these operators read as: 



A* =a - {ij + Gkj] 



^—ik-. 



(2) 



where a is a vector of Dirac matrices and G is an ar- 
bitrary gauge parameter. In the calculations below, fol- 
lowing Grant [35|, we employ two different gauges that 
are known to lead to well known non-relativistic opera- 
tors. First, we use the so-called Coulomb gauge, when G 
—0, which corresponds to the velocity form of electron- 
photon interaction operator in the non-relativistic limit. 
As the second choice we adopt G ~ ^ {L + 1)/L in or- 
der to obtain Babushkin gauge which reduces, for the 
particular case of L=l, to the dipole length form of the 
transition operator. 

In Eq. ll]), \i) = \niKi^i) and |/) = \nfKf^f) denote 
solutions of the Dirac's equation for the initial and final 
ionic states while Ei = E^Ki and Ej = Enj^j are the 
corresponding one-particle energies. Because of energy 
conservation, Ei and Ef are related to the energies Wi^2 = 
cki^2 of the emitted photons by: 



Ei - Ef = ui + u;2 ■ 



(3) 



^From this relation, it is convenient to define the so-called 
energy sharing parameter y ~ u!i/{uji + UJ2), i.e., the 
fraction of the energy which is carried away by the "first" 
photon. 

As usual in atomic physics, the second-order transi- 
tion amplitudes in Eq. ([1]) and, hence, the two-photon 
transitions rates can be further simplified by applying 
the techniques of Racah's algebra if all the operators are 
presented in terms of spherical tensors and if the (stan- 
dard) radial-angular representation of Dirac's wavefunc- 
tions are employed. For the interaction of electron with 
electromagnetic field, the spherical tensor components 
are obtained from the multipole expansion of the oper- 
ator A* (see Refs. [ISj, 19, 39] for further details). By 
using such an expansion, we are able to re-write Eq. ([T]) 
as a sum of partial multipole rates 



dw 
dwi 



E 

eiLie2i2 



dWStLi02L2 

dwi 



(4) 
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which describe the emission of two photons of electric 
{Qj — E) or/and magnetic (9j = M) type carrying away 
the angular momenta Li and L2. For the decay of unpo- 
larized ionic state \niKi), in which the emission angles as 
well as polarization of both photons remain unobserved, 
these partial multipole rates are given by [l8l |: 



L1Q2L2 



duji 



(27r)3c3 ^ ^ 



|5> (1,2)1 



\S=-{2,l)f + 2Y,d{ju,f,)S^"{2,l)S='^{l,2) 



(5) 



where the angular coefficient d{j^,jl) is defined by the 
phase factor and 6j Wigner symbol: 



1 



L2 



(6) 



and the radial integral part is expressed in terms of the 
reduced matrix elements of the multipole (electric and 
magnetic) field operators: 



^^""(1,2 



E 



^ Ae, * 



, Ae, 



El, — Ei -\- UJ2 



■(7) 



The summation over Ae^ in Eq. ([5]) is restricted to Ae^ = 
±1 for the electric {Qj — E) and Ae^ ~ for the magnetic 
{Qj — M) photon transitions. 

Until now, we have discussed the general expressions 
for the two-photon transition rates which are differential 
in energy uji of one of the photons. By performing an in- 
tegration over this energy one may easily obtain the total 
rate that is directly related to the lifetime of a particular 
excited state against the two-photon decay. As it follows 
from Eq. such a total rate can be represented as a 
sum of its multipole components: 



eiLie2L2 

- E 

BiLi02L2 Q 
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dW^ei£ie2-L2 
duji 



dwi , (8) 



where Ut = Ei — Ef is the transition energy. 

As seen from Eqs. (HJ-®, any analysis of the differ- 
ential as well as total two-photon decay rates can be 
traced back to the (reduced) matrix elements that de- 
scribe the interaction of an electron with the (multipole) 
radiation field. Since the relativistic form of these ma- 
trix elements is applied very frequently in studying the 
various atomic processes, we shall not discuss here their 
evaluation and just refer the reader for all details to refer- 
ences [3, [1^ [3^ . Instead, in the next section we will fo- 
cus on the summation over the intermediate states Iki/K^) 
which appears in the second-order transition amplitudes 
(see Eq. ©). 
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FIG. 1: (Color online) Energy-differential transition rates for 
the 2E1 2si/2 lsi/2 two-photon decay of hydrogen and 
hydrogen-like ions. Relativistic calculations have been car- 
ried out by performing intermediate-state summation over 
complete Dirac's spectrum (solid line) as well as by restricting 
this summation to the positive- (dashed line) and negative- 
energy (dotted line) states only. 



B. Summation over the intermediate states 

The summation over the intermediate states in Eq. ([7]) 
runs over the complete one-particle spectrum \nuK.„), in- 
cluding a summation over the discrete part of the spec- 
trum as well as an integration over the positive and 
negative-energy continuum. In practice, of course, per- 
forming such an infinite-state summation is a rather de- 
manding task. A number of methods have been de- 
veloped over the last decades in order to evaluate the 
second-order transition amplitudes consistently. Apart 
from the Green's function approach [l^, [s^ which — 
in case of a purely Coulomb potential — allows for 
the analytical computation of Eq. the discrete- 

basis-set summation is widely used nowadays in two- 
photon studies 19]. A great advantage of the latter 
method is that it allows to separate the contributions 
from the positive- and negative-energy solutions in the 
intermediate-state summation. Since the effects that 
arise from the negative-energy spectrum are in the focus 
of the present study, we apply for the calculations below 
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the finite (discrete) basis solutions constructed from the 
B-sphne sets. 

Although the B-spline basis set approach has been dis- 
cussed in detail elsewhere [l^, [13, [11], here we briefly 
recall its main features. In this way, we shall consider 
the ion (or atom) under consideration to be enclosed in a 
finite cavity with a radius R large enough to get a good 
approximation of the wavefunctions with some suitable 
set of boundary conditions, which allows for discretiza- 
tion of the continua. Wavefunctions that describe the 
quantum states |i/) = \n^Ki,) of such a "particle in box" 
system can be expanded in terms of basis set functions 
(j)l{r) with i — 1, .., 2N which, in turn, are found as so- 
lutions of the Dirac-Fock equation, 



V{r) 



d_ 

dr 



R 

y(r) 



(9) 



where = E], — mc^ and V{r) is a Coulomb potential 
of a uniformly charged finite-size nucleus. Due to com- 
putational reasons, each of 4'\,{r) function is expressed 
as a linear combination of B-splines as it was originally 
proposed in Ref. [13 by Johnson and co-workers. 

For each quantum state \v) the set of basis functions 
<j)\,{r) spans both positive and negative energy solutions. 
Solutions labeled by i = 1,-.,^^ describe the negative 
continuum with e\, < —2mc^ while solutions labeled by 
i = N + \, ..^ 2N correspond to the first few states of the 
bound-state spectrum as well as to positive continuum 
with el > 0. Thus, by selecting the proper sub-set of 
basis functions 4>l,{r) we may explore the role of negative 
continuum in computing of the properties of two-photon 
emission from hydrogen-like ions. 



C. Semi— relativistic approximation 
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FIG. 2: (Color online) Energy-differential decay rates for the 
(sum of the) 2M1, 2E2 and E2M1 2si/2 Isj/j multipole 
two-photon transitions in hydrogen and hydrogen-like ions. 
Relativistic calculations have been carried out by perform- 
ing intermediate-state summation over complete Dirac's spec- 
trum (solid line) as well as by restricting this summation to 
the positive- (dashed line) and negative-energy (dotted line) 
states only. Results of relativistic calculations are compared 
also with the semi-relativistic prediction (dot-dashed line) as 
given by Eq. (fT4|l . 



Based on the relativistic theory, the expressions ob- 
tained in the previous section allow to study the in- 
fiuence of the Dirac's negative continuum on the prop- 
erties of two-photon emission from hydrogen-like ions 
with nuclear charge in the whole range 1 < Z < 92. 
For the low-Z ions, moreover, it is also useful to esti- 
mate the negative-energy contributions within the semi- 
relativistic approach as proposed in the work by Lab- 
zowsky and co-workers [34| . To perform such a semi- 
relativistic analysis let us start from Eq. ^ in which 
we retain the sum only over the negative-energy con- 
tinuum states. Since the total energy of these states is 
~ —{Ti, -t- mc^), the corresponding energy denomi- 
nator of the second-order transition amplitude can be 
written as E^ — Ei + ujj « — 2mc^ which leads to the 
following expression for the differential decay rate: 



dwi 
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(10) 



For the further simplification of this expression we shall 
make use of the multipole expansion of the electron- 
photon interaction operators For the sake of sim- 
plicity, we restrict this semi-relativistic analysis to the 
case of Coulomb gauge {G = 0) in which operator A* 
can be written as: 



a. ■ {I -ik-r + 1/2 {-ik ■ r) + ...) . 



(11) 



if one expand the photon exponential exp(ifc • r) into the 
Taylor series. 

In contrast to the "standard" spherical tensor expan- 
sion [H, [1§| , the series pT|) usually does not allow one 
to make a clear distinction between the different mul- 
tipole components of the electromagnetic field. For in- 
stance, while the first term in Eq. (|lip describes — within 
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the non-relativistic limit — electric dipole (El) transition, 
the term {—ik ■ r) gives rise both, to magnetic dipole 
(Ml) and electric quadrupole (E2) channels. Such an 
approximation, however, is well justified for our (semi- 
relativistic) analysis which just aims to estimate the role 
of negative continuum states in the different [groups of) 
multipole two-photon transitions in light hydrogen-like 
ions. In particular, by adopting A* = —a ■ ej (ik ■ r) for 
both operators in Eq. Ijl0[) we may find the contribution 
from the negative spectrum to the 2M1, 2E2 and E2M1 



2s 



1/2 



lsi/2 transition probabilities: 
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Here, summation over the intermediate states \v) is re- 
stricted by the negative-energy solutions of the Dirac 
equation for the electron in the field of nucleus. In the 
non-relativistic limits these states form a complete set of 
solutions of the Schrodinger equation for the particle in 
a repulsive Coulomb filed |,4Q] . By employing a closure 
relation for such a set we re-write Eq. in the form: 



(-) 

M1.E2 



UJ1UJ2 
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(27r)3c3 (toc^ 



X Keiea) (/ |(fci • r){k2 ■ r)| ^)|' dn^i , (13) 

where \i) and |/) denote now the solutions of the 
Schrodinger equation for the initial and final ionic states, 
respectively. For the particular case of 2si/2 lsi/2 
two-photon transition, i.e., when \i) = \2s) and |/) — 
I Is), this expression finally reads: 
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if one performs an integration over the photon emission 
angles as well as a summation over the polarization states 
(see Ref. [sl] for further details). 

Eq. mi) provides the differential rate for the 2M1, 2E2 
and E2M1 two-photon transitions as obtained within the 
non-relativistic framework and by restricting the summa- 
tion over the intermediate spectrum \i>) to the negative 
energy states only. Being valid for \ow-Z ions, this ex- 
pression may also help us to analyze the negative-energy 
contribution to the total decay rate. 
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FIG. 3: (Color online) Energy-differential decay rates for 
the (sum of the) ElMl and E1E2 2pi/2 lsi/2 multipole 
two-photon transitions in hydrogen and hydrogen-like ions. 
Relativistic calculations have been carried out by perform- 
ing intermediate-state summation over complete Dirac's spec- 
trum (solid line) as well as by restricting this summation to 
the positive- (dashed line) and negative-energy (dotted line) 
states only. Results of relativistic calculations are compared 
also with the semi-relativistic prediction (dot-dashed line) as 
given by Eq. (|16p . 



where the integration over the photon energy uji is per- 
formed. 

Apart from the 2M1, 2E2 and E2M1 2si/2 lsi/2 
two-photon transitions, Eqs. (jlOp and (fTTj) may also be 
employed to study other decay channels. For example, 
the negative energy contributions to the differential as 



well as total rates for the ElMl and E1E2 2pi/2 
decay read as: 



Is 
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and 



"E1,M1,E2 



5.822 10"''^ (aZ)^ 



(17) 



respectively [33] • Together with Eqs. mi) and we 
shall later use these non-relativistic predictions in order 
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to check the validity of our numerical calculations in low- 
Z domain. 



III. RESULTS AND DISCUSSION 

Having discussed the theoretical background for the 
two-photon studies, we are prepared now to analyze the 
influence of the Dirac's negative continuum on the to- 
tal as well as energy-differential decay rates. We shall 
start such an analysis from the 2si/2 — > lsi/2 transi- 
tion, which is well established both in theory [l^j Ej [13 
and in experiment. For all hydrogen-like ions this tran- 
sition is dominated by the 2E1 decay channel while all 
the higher multipoles contribute by less than 0.5% to 
the decay probability. The energy-differential decay rate 
given by Eq. ([5]) for the emission of two electric dipole 
photons is displayed in Fig. [T] for the decay of neutral 
hydrogen (H) as well as hydrogen-like xenon Xe^'^"'" and 
uranium U^'^+ ions. For these ions, relativistic second- 
order calculations have been done within the Coulomb 
gauge and by performing intermediate-state summation 
over the complete Dirac's spectrum (solid line) as well 
as over the positive- (dashed line) and negative-energy 
(dotted line) solutions only. As seen from the figure, the 
negative-energy contribution to the energy-differential 
decay rate is negligible for low-Z ions but becomes rather 
pronounced as the nuclear charge Z is increased. For the 
2E1 decay of hydrogen-like uranium, for example, exclu- 
sion of the negative solutions from the intermediate-state 
summation in Eq. ([7]) leads to about 20 % reduction of 
the decay rate when compared with the "exact" result. 

While for the leading, 2E1 2si/2 — > lsi/2 transition 
the negative continuum effects arise only for rather heavy 
ions, they might strongly affect properties of the higher 
multipole decay channels in low-Z domain. In Fig. [21 
for example, we display the energy distributions of pho- 
tons emitted in 2M1 and 2E2 transitions. As seen from 
the upper panel of the figure corresponding to the decay 
of neutral hydrogen, negative energy part of the Dirac's 
spectrum gives the dominant contribution to the (sum of 
the) differential rates for these decay channels. With the 
increasing nuclear charge Z , the role of positive energy 
solutions also becomes more pronounced. However, these 
solutions allow one to describe reasonably well the differ- 
ential rates ([5]) only if one of the photons is much more 
energetic than the second one, i.e., when either y < 0.1 
or y > 0.9. For a nearly equal energy sharing {y « 0.5), 
in contrast, accurate relativistic calculations of the 2M1 
and 2E2 rates obviously require summation over both, 
the negative and the positive energy states. 

Apart from the results of relativistic calculations, 
we also display in Fig. [2] the (sum of the) negative- 
energy contributions to the 2M1, 2E2, M1E2 and E2M1 
2si/2 — > lsi/2 transition probabilities as obtained within 
the semi-relativistic approach discussed in Section IH CI 
As expected, for low-Z ions both the relativistic (dot- 
ted line) and semi-relativistic (dot-dashed line) results 



basically coincide and are well described by Eq. (|14p . 
As the nuclear charge Z is increased, however, semi- 
relativistic treatment leads to a slight underestimation of 
the negative-energy contribution to the two-photon (dif- 
ferential) transition probabilities. For the 2si/2 — *■ lsi/2 
decay of hydrogen-like uranium ion, for example, results 
obtained from Eq. (fH)) is about 30 % smaller than the 
corresponding relativistic predictions. 

Up to now, we have been considering the 2si/2 lsi/2 
two-photon decay of the hydrogen-like ions. Apart 
from this — experimentally well studied — transition, re- 
cent theoretical interest has been focused also on the 
2pi/2 l'Si/2 two-photon decay [s^l- Although such 
a channel is rather weak comparing to the leading one- 
photon El transition, its detailed investigation is highly 
required for future experiments on the parity violation 
in simple atomic systems [4lj. A number of calcula- 
tions [33, [121 have been performed, therefore, for the 
transition probabilities of the dominant ElMl and E1E2 
multipole components. In order to discuss the role of 
Dirac's negative continuum in these calculations, we dis- 
play in Fig. [3] the energy-differential rate for the sum 
of the ElMl and E1E2 2pi/2 lsi/2 two-photon tran- 
sitions. Again, the calculations have been carried out 
within the Coulomb gauge for the electron-photon cou- 
pling and for three nuclear charges Z = 1, 54 and 92. As 
seen from the figure, negative-energy summation in the 
second-order transition amplitude ^ is of great impor- 
tance for accurate evaluation of 2pi/2 lsi/2 transition 
probabilities both for low-Z and high-Z ions. That is, 
restriction of the intermediate-state summation to posi- 
tive part of Dirac's spectrum results in an overestimation 
of the ElMl and E1E2 differential in energy decay rates 
by factors of about 2 and 2.5 for the neutral hydrogen 
and hydrogen-like uranium, respectively. 

Similarly to the 2si/2 — *■ lsi/2 multipole transitions, 
we make use of semi-relativistic formulae from Section 
III CI to cross-check our relativistic computations for the 
negative-energy contribution to the ElMl and E1E2 
2pi/2 ~^ lsi/2 decay rates in low-Z domain. Again, 
while for neutral hydrogen both, semi-relativistic (|16p 
and relativistic approximations produce virtually identi- 
cal results, they start to differ as the nuclear charge Z is 
increased. 

So far we have discussed the energy-differential de- 
cay rates both for 2si/2 lsi/2 and 2pi/2 lsi/2 
two-photon transitions. Integration of these rates over 
the energy of one of the photons (see Eq. ([5])) will yield 
the total decay rates. In Table I we display the to- 
tal decay rates for the various multipole channels of 
2si/2 lsi/2 two-photon decay. In contrast to the 
photon energy distributions from above, here relativis- 
tic calculations have been performed in Coulomb (ve- 
locity) as well as Babushkin (length) gauges. In both 
gauges, negative-energy contribution to the (total) prob- 
ability of the leading 2E1 transition is about eight orders 
of magnitude smaller than positive-energy term if de- 
cay of low-Z ions is considered but is significantly in- 
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TABLE I: Total rates (in s~^) for the several multipole combinations of 2si/2 —* lsi/2 two-photon decay. Relativistic calcula- 
tions have been performed within the velocity and length gauges and by carrying out intermediate-state summation over the 
complete Dirac's spectrum (Wt) as well as over the positive- (W+) and negative-energy (W_) solutions only. 



creased for higher nuclear charges. For the hydrogen- 
hke uranium, for example, the total 2E1 decay rate is 
enhanced from 2.9041x10^^ s~^ in the velocity gauge 
and 2.3939x10^^ s~^ in the length gauge to the — gauge 
independent — "exact" value of 3.8256x 10^^ s~^ if, apart 
from the positive-energy states, the Dirac's states with 
negative energy are taken into account in transition am- 
phtude Eq. ([7]). These results clearly indicate the impor- 
tance of the negative-state summation for the accurate 
evaluation of 2E1 2si/2 — > lsi/2 total rates in both, ve- 
locity and length gauges. It worth mentioning, however, 
that while for velocity gauge our findings are in perfect 
agreement with results reported in Ref. some dis- 
crepancy was found for calculations performed in length 
gauge for which Labzowsky and co-workers have argued 
that the contribution from the Dirac's negative contin- 
uum is negligible even for heaviest ions. The reason for 
this discrepancy is not apparent for the moment and, 
hence, further investigations are highly required. 

In Table I, besides the leading 2E1 decay channel, 
we present the results of relativistic calculations for the 
higher multipole contributions to the 2sx/2 ~* l'Si/2 two- 
photon transition. The influence of Dirac's negative con- 
tinuum is obviously different for various multipole com- 
binations. While, for example, the negative-energy con- 
tribution to the intermediate-state summation in low-Z 
domain is negligible for the E1M2 decay it becomes of 
paramount importance for the 2E2 and 2M1 decay chan- 
nels; an effect that has been already discussed for the case 
of the energy-differential decay rates (see upper panel of 
Fig. [2]). Moreover, 2sx/2 ~* lsi/2 transition with emis- 
sion of two magnetic dipole (2M1) photons in light ions 
seems to happen almost exclusively via the negative en- 
ergy (virtual) intermediate states. The total decay rate 
for this transition together with the negative-energy con- 
tribution to the probability of the 2E2 channel (evaluated 



in Coulomb gauge) gives in atomic units: 

W2M1 + W2£;2 = 1-248 10"^ {aZ^^ , (18) 

which is in perfect agreement with the semi-relativistic 
formula 

As mentioned above for the computation of the pho- 
ton energy distributions in \ow-Z domain, negative- 
energy contribution to the intermediate-state summation 
is rather pronounced not only for the higher multipole 
terms of 2si/2 lsi/2 decay but also for the leading 
ElMl and E1M2 (two-photon) channels of 2pi/2 — > lsi/2 
transition. Our relativistic calculations displayed in Ta- 
ble II indicate that one should account for negative- 
continuum summation also for an accurate evaluation of 
the total decay rates for these two decay channels. For 
the decay of light elements, sizable contribution from the 
negative-continuum intermediate states arises both in 
length and velocity gauges. Again, these results partially 
question the predictions by Labzowsky and co-workers 
[34| who claimed a minor role of negative energy terms 
for ElMl and E1M2 calculations in length gauge. For the 
velocity gauge, in contrast, our relativistic calculations: 

^Eiiii + ^EiE2 = 5-822 10-^ (aZ)« , (19) 

are in a good agreement both, with the senii-relativistic 
prediction ([TT]) and data presented in Ref. (S^ . 

IV. SUMMARY AND OUTLOOK 

In conclusion, the two-photon decay of hydrogen-like 
ions has been re-investigated within the framework of 
second-order perturbation theory, based on Dirac's rel- 
ativistic equation. Special attention has been paid to 
the summation over the intermediate ionic states which 
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TABLE II: Total rates (in s ^) for the several multipole combinations of 2pi/2 — > lsi/2 two-photon decay. 



occurs in such a framework and runs over complete one- 
particle spectrum, including a summation over discrete 
(bound) states as well as the integration over the posi- 
tive and negative continua. In particular, we discussed 
the role of the negative energy continuum in an accurate 
evaluation of the second-order transition amplitudes and, 
hence, the energy-differential as well as total decay rates. 
Detailed calculations of these rates have been presented 
for the 2si/2 ~* lsi/2 and 2pi/2 lsi/2 two-photon 
transitions in neutral hydrogen as well as hydrogen-like 
xenon and uranium ions. As seen from the results ob- 
tained, both the total decay probabilities and the en- 
ergy distributions of the simultaneously emitted photons 
can be strongly affected by the negative-state summa- 
tion not only for heavy ions but also for low-Z domain. 
We demonstrate, however, that the role of Dirac's nega- 
tive continuum becomes most pronounced for the higher 
(non-dipole) terms in the expansion of the electron- 
photon interaction; similar effect has been recently re- 
ported for the theoretical description of hydrogen-like 
systems exposed to intense electromagnetic pulses |30l] . 

In the present work, we have restricted our discussion 
of the negative energy contribution to the second-order 
calculations of the total and energy-differential decay 
rates. Even stronger effects due to the Dirac's negative 



continuum can be expected, however, for the angular and 
polarization correlations between emitted photons. The- 
oretical investigation of these correlations which requires 
also detailed analysis of interference terms between the 
various (two-photon) multipole combinations is currently 
underway and will be published soon. 
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